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Abstract. This paper presents an algorithm that achieves hyper-arc 
consistency for the soft alldifferent constraint. To this end, we prove 
and exploit the equivalence with a minimum-cost flow problem. Consis- 
tency of the constraint can be checked in 0(nm) time, and hyper-arc 
consistency is achieved in 0(m) time, where n is the number of variables 
involved and m is the sum of the cardinalities of the domains. It improves 
a previous method that did not ensure hyper-arc consistency. 

1 Introduction 

If a constraint satisfaction problem (CSP) is over-constrained, i.e. has no so- 
lution satisfying all constraints, it is natural to allow certain constraints, the 
soft constraints, to be violated and search for solutions that violate as few soft 
constraints as possible. Constraints that are not decided to be soft are hard 
constraints, and should always be satisfied. 

Several methods have been proposed to handle over-constrained CSPs, see for 
instance [6,2,4]. In this paper, we follow the scheme proposed by Regin, Petit, 
Bessiere and Pugct [11], that is particularly useful for non-binary constraints. 
The idea is as follows. A cost function is assigned to each soft constraint, measur- 
ing the violation. Then the soft CSP is transformed into a constraint optimiza- 
tion problem (COP), where all constraints are hard, and the (weighted) sum of 
cost functions is minimized. This approach allows one to use specialized filtering 
algorithms for soft constraints, as shown by Petit, Regin and Bessiere [7]. 

For the soft alldifferent constraint, an algorithm is presented in [7] that re- 
moves inconsistent values in 0(m 2 n^/n) time, where n is the number of variables 
and m the sum of the cardinalities of their domains. However, that algorithm 
does not ensure hyper- arc consistency. In this paper, we propose an algorithm 
that does ensure hyper-arc consistency and runs in O(nm) time. In principle, we 
consider the soft alldifferent constraint as a minimum-cost flow problem in 
a particular graph. Checking the consistency can then be done in 0{nm) time. 
Thereafter, domain values are checked for consistency by an efficient shortest 
path computation, which takes in total O(m) time. 

The outline of the paper is as follows. Section 2 presents definitions related 
to constraint satisfaction problems. Section 3 shows a graph-theoretic analysis 



of the soft alldif f erent constraint, using flow theory. In Section 4 the filtering 
algorithm is presented. We conclude with a discussion in Section 5. 

2 Background 

We assume familiarity with the basic concepts of constraint programming. For 
a thorough explanation of constraint programming, see [1]. 

A constraint satisfaction problem (CSP) consists of a finite set of variables 
V = {«!,..., v r } with finite domains T> = {Di,...,D r } such that Vi E Di 
for all i, together with a finite set of constraints C, each on a subset of V. A 
constraint C E C is defined as a subset of the Cartesian product of the domains 
of the variables that are in C. A tuple [d\, . . . , d r ) E D\ x • • • x D r is a solution 
to a CSP if for every constraint C E C on the variables v i± , . . . , Vi k we have 
(diu . . . ,d ik ) e C. A constraint optimization problem (COP) is a CSP together 
with an objective function to be optimized. A solution to a COP is a solution 
to the corresponding CSP, that has an optimal objective function value. 

Definition 1 (Hyper-arc consistency). A constraint C on the variables x\, 
. . . ,x k is called hyper-arc consistent if for each variable Xi and value di E Di, 
there exist values d\, . . . , dj_i, d i+1 , . . . , d k in D\, . . . , D i+ i, . . . , Df~, such 

that (di, . . . , dk) <G C. 

Definition 2 (Consistent CSP). A CSP is hyper-arc consistent if all its con- 
straints are hyper- arc consistent. A CSP is inconsistent if it has no solution. 
Similarly for a COP. 

Definition 3 (Pairwise difference). Let X\,...,x n be variables with respec- 
tive finite domains Di , . . . , D n . Then 

alldif ferent(xi, . . . ,x n ) = {(di, . . . ,d n ) \ di E Di,dj ^ du for j ^ k}. 

In [7], two different measures of violation for a soft constraint are presented. 
The first is the minimum number of variables that need to change their value 
in order to satisfy the constraint. For this measure, applied to the alldiff- 
erent constraint, [7] also contains a hyper-arc consistency algorithm. The second 
measure is the number of violated constraints in the binary decomposition of the 
constraint, if this decomposition exists. For the alldif f erent constraint, such a 
decomposition does exist, namely Xi ^ Xj for i E {1, . . . , n— l},j € {i+l, . . . ,n}. 
We follow this second, more refined, measure, and present it in terms of the 
soft alldif f erent constraint. For alldif ferent(xi, x n ), let the cost of 
violation be defined as 

violation(xi, ...,x n ) = \ x { = Xj, for i < j}\ . (1) 

Definition 4 (Soft pairwise difference). Let . . . , x n , z be variables with 
respective finite domains D\, . . . , D ni D z . Then 

sof t_alldif f erent(xi, . . . , x n , z) = 

{{d\, . . . ,d n ,d) | di e Di, d E D z , violation(di, . . . , d n ) < d}. 



The variable z in Definition 4 will serve as a so-called cost variable, which will 
be minimized during the solution process. This means that admissible tuples 
in Definition 4 are those instantiations of variables, such that the number of 
violated dis-equality constraints di ^ dj is not more than that of the currently 
best found solution, represented by maxD z . At the same time, min_D z should 
not be less than the currently lowest possible value of violation^!, . . . , x n ). 

An over-constrained CSP with an alldif f erent constraint is transformed 
into a COP by introducing z, replacing alldif f erent with sof t_alldif f erent, 
and minimizing z. This is illustrated in the following example. 

Example 1. Consider the following over-constrained CSP 

x\ £ {a, b}, x 2 £ {a, 6}, x 3 g {a, b}, x 4 g {b, c}, 
alldif f erent(xi, x 2 , Xs, x±). 

We transform this CSP into 

«e{0 6}, 

x\ g {a, 6}, x 2 £ {a, b}, x 3 £ {a, b}, x 4 £ {b, c}, 
sof t_alldif f erent(a;i, x 2 , x 3 ,x 4 , z), 
minimize z. 

This COP is not hyper-arc consistent, as there is no support for z < 1. If we 
remove from D z , the COP is hyper-arc consistent, because there are at most 6 
simultaneously violated dis-equalities. Suppose now that during the search for a 
solution, we have found the tuple (xi, x 2 , x$, x 4 , z) = (a, a, 6, c, 1), that has one 
violated dis-equality. Then z £ {1} in the remaining search. As the assignment 
x 4 = b always leads to a solution with z > 2, b can be removed from D 4 . The 
resulting COP is hyper-arc consistent again. 

One should take into account that a simplified CSP is considered in Exam- 
ple 1. In general, a CSP can consist of many more constraints, and also more 
cost-variables that together with z form an objective function to be minimized. 

Throughout this paper, let m = J2ie{i «} 1-^*1 ^ or variables x\, . . . , x n . 

3 Graph-Theoretic Analysis 

A directed graph is a pair G = (V, A) where V is a finite set of vertices V and 
A is a family 1 of ordered pairs from V, called arcs. A pair occurring more than 
once in A is called a multiple arc. For v £ V, let S m (v) and <5 out (v) denote the 
family of arcs entering and leaving v respectively. 

A (directed) walk in G is a sequence P = vo, cii, v\, . . . , a,k, Vk where fc > 0, 
Vo, Vi, . . . , Vk £ V, ai, . . . , flfc g A and a, = u,) for i = 1, . . . , k. If there 

is no confusion, P may be denoted as P — v 0} vi, . . . , v k . A (directed) walk is 
called a (directed) path if v 0} . . . ,v k are distinct. A closed (directed) walk, i.e. 
vo = Vk, is called a (directed) circuit if v\, . . . , v k are distinct. 

1 A family is a set in which elements may occur more than once. 



3.1 Minimum-cost flow problem 



First, we introduce the concept of a flow, following Schrijver [12, pp. 148-150]. 

Let G = (V, A) be a directed graph and let s, t e V. A function / : A — > R is 
called a flow from s to t, or an s — t flow, if 

{i) /(a) > for each oei, . . 

(ii) f(5 out (v)) = f(S in (v)) for each v e V \ {s, i}, 

where /(S) = X) a es-^( a ) ^ or au ^ — Property (2)(m) ensures flow conser- 
vation, i.e. for a vertex v ^ s,t, the amount of flow entering v is equal to the 
amount of flow leaving v. 

The value of an s — t flow / is defined as 

value(/) = f(6™\s)) - f(S in (s)). 

In other words, the value of a flow is the net amount of flow leaving s, which 
can be shown to be equal to the net amount of flow entering t. 

When we study flows we typically endow capacity constraints, via a "capac- 
ity" function c : A — > R + . We say that a flow / is under c if /(a) < c(a) for each 
a e A. A feasible flow is a flow under c. 

We also assign costs to flows via a "cost" function w : A — > K + . Doing so the 
cost of a flow / is defined as 

cost(/) = w(a)f(a). 

A minimum-cost flow is an s — t flow under c of maximum value and minimum 
cost. The minimum-cost flow problem is the problem of finding such a minimum- 
cost flow. 

A minimum-cost flow can be computed using an algorithm originally due 
to Ford and Fulkerson [5] (we follow the description given by Schrijver [12, pp. 
183-185]). It consists of successively finding shortest (with respect to the cost 
function) s—t paths in the so-called residual graph, while maintaining an optimal 
flow. 

Define the residual graph Gf = (V, Af) of / (with respect to c), where 
A f = {a | a e A, f(a) < c(o)} U {aT 1 | a e A, f(a) > 0}. 
Here a -1 = (v, u) if a = (u, v). We extend to to A^ 1 = {a -1 | a e A} by defining 

t#(a _1 ) = —w(a) 

for each a £ A. 

Any directed path P in Gf gives an undirected path in G = (V, A). We define 
X p eR A by 

!1 if P traverses a, 
— 1 if P traverses a -1 , 
if P traverses neither a nor a -1 , 



Algorithm 1 Minimum-cost s — t flow 



set / = 

while termination criterion not satisfied do 
compute minimum-cost s — t path P in Gf 
if no s — t path in Gf exists then 

terminate 
else 

set e maximal, such that < / + ex P < c 
reset f = f + e\ P 
end if 
end while 



for a £ A. Define the cost of a path P as cost(P) = J2aeP w ( a )- 

Call a feasible flow extreme when it has minimum cost among all feasible 
flows with the same value. Then the following holds (cf. [12, Theorem 12.3 and 
12.4]). Let denote the all-zero vector of appropriate size. 

Theorem 1. A flow f is extreme if and only if each directed circuit of Gf has 
nonnegative cost. 

Theorem 2. Let f be an extreme flow in G = (V, A). Let P be a minimum-cost 
s — t path in Gf, for some s,t € V, and let e > be such that f = f + ex P 
satisfies < /' < c. Then f is an extreme flow again. 

In fact, for /, P, e and /' in Theorem 2 holds 

value(/') = value(/) + e, 
cost(/') = cost(/) + e • cost(P). 

This means that we can find a minimum-cost s — t flow in G by successively 
computing minimum-cost s — t paths in Gf. Along such a path we increase the 
amount of flow to the maximum possible value s. By Theorem 2, the last flow 
(of maximum value) we obtain must be extreme, and hence optimal. This is 
presented as Algorithm 1 . Note that the cost of minimum-cost s — t paths in G / 
is bounded, because there are no directed circuits of negative cost in Gf. For 
rational capacities, Algorithm 1 terminates with a feasible s — t flow of maximum 
value and minimum cost. Although faster algorithms exist for general minimum- 
cost flow problems, Algorithm 1 suffices when applied to our problem. This is 
because in our particular graph Algorithm 1 is faster than the algorithms for 
general minimum-cost flow problems. 

3.2 From sof t_alldif f erent to minimum-cost flow 

We transform the problem of finding a solution to the sof t_alldiff erent con- 
straint into a minimum-cost flow problem. 

Construct the directed graph G = (V, A) with 



V = {a, t}UXUD x 



Fig. 1. Graph G for the sof t_alldif f erent constraint of Example 1. For each arc a, 
(c(a), w(a)) is given. Bold arcs indicate an optimal s — t flow with cost 1. 



and 

A = A x U A s U A t 

where 

X = {xi, . . .,x n }, 

D X = Ui£{l,..,n}A, 

and 

A x = {(xud) \deD,}, 
A s = {{s,xi) | i G {l,...,n}}, 
A t = {(d,t) | deDi,i e {l,...,n}}. 

Note that A t contains parallel arcs if two or more variables share a domain value. 
If there are k parallel arcs (d, t) between some d e Dx and t, we distinguish them 
by numbering the arcs as (d, t)o, (d, t)i, . . . , (d, t)k-i in a fixed but arbitrary way. 

To each arc a e A, we assign a capacity c(a) = 1 and a cost w(a). If a e 
A s U Ax, then ty(a) = 0. If a e A t , so a = for some d G -Dx and integer 

i, the value of w(a) = i. 

In Figure 1, the graph G for the sof t_alldiff erent constraint in Example 1 
is depicted. For each arc a, (c(a),w(a)) is given. 

Theorem 3. An integer flow f that is a solution to the minimum-cost flow prob- 
lem in G corresponds to an instantiation of variables X\,...,x n in soft_all- 
diff erent (ari, . . . ,x n ,z), minimizing violation(xi, . . . ,x n ). 

Proof. For an integer flow / in G, f(a) = 1 if arc a is used, and f(a) = oth- 
erwise. An arc a = (xi,d) € Ax with /(a) = 1 corresponds to the instantiation 
Xi = d. By construction, every solution / to the minimum-cost flow problem in 
G has value (/) = n. Thus a solution corresponds to assigning a value to each 
variable x%, i € {1, . . . , n}. 

The cost function io(aj) = i for k parallel arcs oo, ...,afc_i € A t corre- 
sponds to counting the number of violations caused by assigning i + 1 vari- 
ables to a particular value. Namely, for these parallel arcs, a minimum-cost 
s — t path in Gf uses the arc with lowest cost first. Using arc ctj (the (i + l)st 



arc) causes a "violation" with the i previously used arcs. Thus, for a feasi- 
ble flow /, which corresponds to an assignment of x\,.. . ,x n , J2 a eA w ( a )f( a ) 
measures exactly violation (x i, . . . ,x n ). Hence, a minimum-cost flow minimizes 
violation(a;i , . . . , x n ). □ 

Consider again the graph G in Figure 1. A bold arc a in G denotes f(a) = 1. 
This particular flow / has value (/) = 4 and cost(/) = f. Indeed, the only 
violation is X\ = a = x 2 - 

Next we describe the behaviour of Algorithm f to compute a minimum-cost 
flow in G. We need to compute a sequence of minimum-cost s — t paths in G/ , 
maintaining extreme intermediate flows. Note that along each minimum-cost s—t 
path in G / we can increase the flow by a maximum of e = 1 . Hence all extreme 
flows in G arc integer. By construction, there are exactly n such paths, each 
containing one arc in A s (in fact, the paths may as well be computed starting 
from the vertices Xi instead of s, using only arcs in Ax and A t ). Further, each 
minimum-cost s — t path contains exactly one arc in A t . Namely, consider a 
minimum-cost path P using multiple arcs in A t . Then P consists of an s — t path 
with one arc in A s , followed by a t— t path. If the t—t path has cost 0, we may omit 
this part, and use only the s — t path with one arc in A s . If the t — t path, which 
is a circuit, has negative cost, it contradicts Theorem 1. Effectively, it means 
that the t — t path could have been used to improve the preceding intermediate 
solution, thus contradicting the extremity of that solution. To conclude, the 
minimum-cost paths we need to compute use exactly one arc in A s and one arc 
in A t . It follows that these paths can be computed in 0(m) time, and the total 
time complexity for finding a maximum flow of minimum cost in G is 0(nm). 
Hence it follows, by Theorem 3, that consistency of the sof t_alldif f erent 
constraint can be checked in 0{nm) time. 



4 The Filtering Algorithm 

The following theorem identifies hyper-arc consistent domain values for the 
sof t_alldiff erent constraint. For an arc a of G, let G a arise from G by en- 
forcing f(a) = 1 for every flow / in G. 

Theorem 4. The constraint sof t_alldif f erent (xi, . . . , x n , z) is hyper-arc con- 
sistent if and only if 

(i) for all all arcs a G Ax a minimum-cost flow of maximum value in G a has 

cost at most maxfl z , 
(ii) all values in D z are not smaller than the cost of a minimum-cost flow of 

maximum value in G. 

Proof. Enforcing f(a) = 1 for arc a — (xi,d) corresponds to assigning Xi — d. 
The result follows from Definition 1 and Theorem 3. Namely, property (i) checks 
consistency for all domain values in D\, . . . , D n . Property (ii) checks consistency 
of the domain values of D z . □ 



Algorithm 2 Naive hyper-arc consistency 
set minimum = oo 
for Xi € X do 
for d 6 ft do 

compute minimum-cost s — t flow / in G a where a — (xi, d) 
if cost(/) > maxft then 

remove d from Di 
end if 

if cost(/) < minimum then 

set minimum = cost(/) 
end if 
end for 
end for 

if min D z < minimum then 

set min D z = minimum 
end if 



Using Theorem 4, we can construct an algorithm that enforces hyper-arc 
consistency for the sof t_alldif f erent constraint, presented as Algorithm 2. 
For all variables ij £ I, the algorithm scans all domain values d £ Di, and 
checks whether there exists a minimum-cost s — t flow in G a , where a = (xi,d), 
of maximum value with cost at most m&xD z . If such a flow does not exist, then, 
by Theorem 4, d is removed from Di. Finally, we remove all values of D z which 
are smaller than the cost of a minimum-cost flow in G. The time complexity of 
Algorithm 2 is 0(m 2 n). 

We can construct a more efficient filtering algorithm, however. It is presented 
as Algorithm 3, and makes use of the following theorem. We follow the notation 
introduced in Section 3.1. 

Theorem 5. Let f be an extreme flow of maximum value in G. Let a = (xi,d) £ 
Ax and P a minimum-cost d — Xi path in Gf. Let f* be an extreme flow of 
maximum value in G a . Then cost(/*) = cost(/) +cost(P). 

Proof. Either f(a) = 1 or /(a) = 0. In case /(a) = 1, /*(a) = 1, P = d,Xi, 
cost(P) = and we are done. In case /(a) = 0, first note that there exists a 
d — Xi path in Gf. Namely, there is exactly one d! € Di for which f((xi,d')) = 1, 
which allows the path d, t, d' , Xj. Let P be a minimum-cost d — x, path in Gf. 
Together with arc (xj, d) P forms a circuit C. The directed circuit C in Gf gives 
an undirected circuit in G a . For all 6 € A, define flow /* in G a as follows: 

f*(b) = } 1 if be C 
[ f(b) else. 

It is easy to check that /* is again a flow of maximum value. 

Because / is extreme, we may assume that P enters and leaves t only once, say 
via arcs b ln and 6 out respectively (where b m = (d,t)). It follows that cost(P) = 



Algorithm 3 More efficient hyper-arc consistency 
compute minimum-cost flow / in G 
if cost(/) > maxZ>2 then 

return Inconsistent 
end if 

if minD z < cost(/) then 

set minD z = cost(/) 
end if 

for a = (xi,d) with f(a) — do 

compute minimum-cost d — Xi path P in Gf 
if cost(/) + cost(P) > maxDz then 

remove d from Di 
end if 

end for 



w(bin) — w(b out )- From Theorem 1 we know that cost(P) > 0. Similarly, 

cost(r) = £ t6Jl /*(&M6) 

= EfceA f( b ) w ( b ) + w(b in ) - w(b out ) 
= cost(/) + cost(P) 

It remains to show that /* is extreme in G a . Suppose not, i.e. there exists 
a flow g in G a with maximum value and cost(g) < cost(/*). As cost(/*) = 
cost(/) + cost(P) and cost(P) > 0, there are two possibilities. The first is that 
cost(g) < cost(/), which is not possible because / is extreme. The second is that 
there exists an Xi — d path P' in Gf with cost(P') < cost(P) which also leads 
to a contradiction because P is a minimum-cost path. Hence /* is extreme. □ 

Algorithm 3 first computes a minimum-cost flow / in G. This takes 0{nm) 
time, as we have seen in Section 3.2. If cost(/) > maxP/ z , we know that the 
sof t_alldif f erent constraint is inconsistent. If this is not the case, we update 
minl? z . Next, we scan all arcs a — (xi, d) for which /(a) = 0. For each of these 
arcs, we compute a minimum-cost d — Xi path P in Gf. By Theorem 5 and 
Theorem 4, we remove d from Di if cost(/) + cost(P) > maxD z . This can be 
done efficiently, as shown by the following theorem. 

Theorem 6. Let sof t_alldif f erent(.Ti, . . . ,x n ,z) be consistent and f an in- 
teger minimum- cost flow in G. Then sof t_alldif f erent(a;i, . . . , x n , z) can be 
made hyper- arc consistent in 0(m) time. 

Proof. The complexity of the filtering algorithm depends on the computation 
of the minimum-cost d — Xi paths for arcs (xi,d). We make use of the fact that 
only arcs a <G At contribute to the cost of such a path. 

Consider the strongly connected components 2 of the graph Gf which is a copy 
of Gf where s and t and all their incident arcs are removed. Let P be a minimum- 
cost d— Xi path P in Gf. If P is equal to d, Xi then f(xi, d) — 1 and cost(P) = 0. 

2 A strongly connected component in a directed graph G = (V, A) is a subset of vertices 
S C V such that there exists a directed u — v path in G for all u,v G S. 



Otherwise, either xi and d are in the same strongly connected component of Gf, 
or not. In case they are in the same strongly connected component, P can avoid 
t in Gf, and cost(P) = 0. In case Xi and d are in different strongly connected 
components of G f , say Xi € S\ and d G S2 , we have 

cost(P) = min w(a) + min w(a). (3) 

a e {(<i',t) I (d',t) e A f , a e {(t.d") | (t,d") G A f , 

d' £ S 2 } d" £ Si or (d",Xi) £ Af} 

Property (3) follows from the fact that P uses exactly one ingoing and one 
outgoing arc for t. 

Arcs a with f(a) = 1 or within a strongly connected component will all use 
a minimum-cost path with cost equal to 0, and will therefore be all consistent if 
cost(/) < maxD z . For all other arcs, we can resort to property (3). For this we 
only have to compute once for each strongly connected component S of Gf the 
minimum-cost arc going from S to t and the minimum-cost arc going from t to 
S (if such arcs exist), which takes in total 0(m) time. The strongly connected 
components of Gf can be computed in 0(n + m) time, following Tarjan [14]. 
Hence the total time complexity of achieving hyper- arc consistency is 0(m), as 
n < m. □ 

The proof of Theorem 6 applies to any constraint whose graph representation 
resembles G and has only costs on arcs from Dx to t. For all such constraints 
that are consistent, hyper-arc consistency can be achieved in O(m) time. Note 
that this is equal to the complexity of achieving hyper-arc consistency on these 
constraints if no costs are involved. 

5 Conclusion and Discussion 

We have presented an algorithm that checks consistency of the sof t_alldif f- 
erent constraint on n variables in 0(nm) time and achieves hyper-arc consis- 
tency in 0(m) time, where m is the sum of the cardinalities of the domains. 
A previous method for removing domain values that are inconsistent with the 
sof t_alldif f erent constraint did not ensure hyper-arc consistency [7]. More- 
over, that method has a time complexity of 0{w?n^fn). Hence our algorithm 
improves on this in terms of quality as well as time complexity. 

The sof t_alldiff erent constraint is related to the standard alldiff erent 
constraint [8] and the minimum weight alldiff erent constraint [3]. The mini- 
mum weight alldiff erent constraint is a particular instance of the global car- 
dinality constraint with costs [9, 10]. For that constraint, hyper-arc consistency 
can be achieved in 0(n(m + d\ogd)) time, where d is the cardinality of the union 
of all domains [9, 10, 13]. It is achieved by finding n shortest paths, each taking 
0(m + dlogd) time to compute. Although our algorithm has a similar flavour, 
the underlying graphs have a different cost structure. We improve the efficiency 
by exploiting the cost structure of our particular graph when computing the 
shortest paths. Our result can be applied to other constraints with a similar 
graph representation and cost structure. 



Acknowledgements 



Many thanks to Bert Gerards for valuable comments. Thanks also go to Se- 
bastian Brand for fruitful discussion. Finally, the constructive remarks of Jean- 
Charles Regin were highly appreciated. 

References 

1. K.R. Apt. Principles of Constraint Programming. Cambridge University Press, 
2003. 

2. S. Bistarelli, U. Montanari, and F. Rossi. Semiring-based Constraint Satisfaction 
and Optimization. Journal of the ACM, 44(2):201-236, 1997. 

3. Y. Caseau and F. Laburthe. Solving Various Weighted Matching Problems with 
Constraints. In G. Smolka, editor, Proceedings of the Third International Confer- 
ence on Principles and Practice of Constraint Programming (CP'97), volume 1330 
of LNCS, pages 17-31. Springer, 1997. 

4. M.C. Cooper and T. Schiex. Arc consistency for soft constraints. Artificial Intel- 
ligence, 2004. To appear. 

5. L.R. Ford and DR. Fulkerson. Constructing maximal dynamic flows from static 
flows. Operations Research, 6:419-433, 1958. 

6. E.C. Freuder and R.J. Wallace. Partial constraint satisfaction. Artificial Intelli- 
gence, 58(l-3):21-70, 1992. 

7. T. Petit, J.-C. Regin, and C. Bessiere. Specific Filtering Algorithms for Over- 
Constrained Problems. In T. Walsh, editor, Proceedings of the Seventh Inter- 
national Conference on Principles and Practice of Constraint Programming (CP 
2001), volume 2239 of LNCS, pages 451-463. Springer, 2001. 

8. J.-C. Regin. A Filtering Algorithm for Constraints of Difference in CSPs. In 
Proceedings of the Twelfth National Conference on Artificial Intelligence (AAAI- 
94), volume 1, pages 362-367, 1994. 

9. J.-C. Regin. Arc Consistency for Global Cardinality Constraints with Costs. In 
J. Jaffar, editor, Proceedings of the Fifth International Conference on Principles 
and Practice of Constraint Programming (CP'99), volume 1713 of LNCS, pages 
390-404. Springer, 1999. 

10. J.-C. Regin. Cost-Based Arc Consistency for Global Cardinality Constraints. Con- 
straints, 7:387-405, 2002. 

11. J.-C. Regin, T. Petit, C. Bessiere, and J.-F. Puget. An Original Constraint Based 
Approach for Solving over Constrained Problems. In R. Dechter, editor, Proceed- 
ings of the Sixth International Conference on Principles and Practice of Constraint 
Programming (CP 2000), volume 1894 of LNCS, pages 543-548. Springer, 2000. 

12. A. Schrijver. Combinatorial Optimization - Polyhedra and Efficiency. Springer, 
2003. 

13. M. Sellmann. An Arc-Consistency Algorithm for the Minimum Weight All Different 
Constraint. In P. Van Hentenryck, editor, Proceedings of the Eighth International 
Conference on Principles and Practice of Constraint Programming (CP 2002), 
volume 2470 of LNCS, pages 744-749. Springer, 2002. 

14. R. Tarjan. Depth-first search and linear graph algorithms. SIAM Journal on 
Computing, 1:146-160, 1972. 



